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ABSTRACT 

We present imaging and spectroscopic observations with HST and VLT of the ring of SN 1987A 
from 1994 to 2014. After an almost exponential increase of the shocked emission from the hotspots 
up to day ~ 8,000 2009), both this and the unshocked emission are now fading. From the radial 

positions of the hotspots we see an acceleration of these up to 500 — 1000 kms -1 , consistent with 
the highest spectroscopic shock velocities from the radiative shocks. In the most recent observations 
(2013 and 2014), we find several new hotspots outside the inner ring, excited by either X-rays from 
the shocks or by direct shock interaction. All of these observations indicate that the interaction with 
the supernova ejecta is now gradually dissolving the hotspots. We predict, based on the observed 
decay, that the inner ring will be destroyed by ~ 2025. 

Subject headings: circumstellar matter — shock waves - supernovae: individual (SN 1987A) 
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1. INTRODUCTION 


Supernova 1987A is unique: we can study in real time 
and with high spatial resolution physical processes that 
we can only infer for other supernovae. Emission from SN 
1987A included a flash of ultraviolet radiation in the first 
hours, thermal diffusion of radioactive energy input for 
weeks, and prom pt rep r oces sing of radioactive power for 
decades (e.g., McCray|l9 93). But now, conversion of the 
supernova’s kinetic energy to heat through sho ck inter¬ 
acti ons dominates the radiative energy budget (Larsson 
et al. 2 011 in the following Lll). 


The ring system around SN 1987A was created by mass 
loss from the pr ogenito r star i n the 20,000 years before 
the explosion (Crotts & Heathcote 1991). The rings 
may have been created as a res ulto f mass ejections from 
a rapidly rotating single star ([Chita et a l. 2008) or a 
merger of two massive stars (Morris & Fodsiadlowski 
20091 Akashi et al.|2015 ). Observations of the interaction 


of the ejecta with the circumstellar medium can provide 
clues to the mass loss history and thereby the nature of 
the progenitor. 

A blast wave, driven by the supern ova ejecta, ini tiall y 
with a velocity of ;> 35, 000 kms -1 (Staveley-Smith et al 
1993), slowed down to ~ 4000 kms -1 after interacting 
with the H II region created by the B3 la progenitor 
star in the swept-up red supergiant wind (Chevalier fc 


]|1995| |Dewey et al.||2012| |Fotter et al.||2014|). 
teraction with the circumstellar ring occurred 


Dwarkadas 

T he first interaction witR tlie circumstellar ring 
in 1995. This event was manifested by the appearance 
of the first hotspot (|Sonneborn et al ||1998| Lawrence 
et al. 2000), where the blast wave entered a relatively 
dense clump of matter and suddenly slowed to velocities 
of a few hundred kms -1 . Additional hotspots erupted in 
dense clumps of gas over the next few years, gra dual ly in¬ 
volving the whole ring . Emission at optical (Groningsson 


volving the whole rin g. Emission at optical (Groningsson 
et al.|2008al), X-ray (Maggi et al. 2012jjHelder et al.|2013p 
and radio ( |Ng et ah ' 2013| ) wavelengths have since then 


undergone a steep rise. 

We have followed the evolution of these shocks from 
the very beginning through regular monitoring with the 
Hubble Space Telescope (HST) and the Very Large Tele¬ 
scope (VLT). In this Letter we report observations of the 
ring system over the last ^ 20 years, and show that the 
ring interaction has entered a new phase involving the 
final destruction of the ring but also interactions with 
material outside of this. 


2. OBSERVATIONS AND ANALYSIS 

In order to make light curves of the ring we use all 
HST imaging observations in the R-band (F675W and 
F625W) and B-band (F439W, F435W and F438W) ob¬ 
tained since 1994. Details of all observations up to De¬ 
cember 2009 are included in Lll. Additional observa¬ 
tions are given in Table [l] All images were processed in 
a standard way in order to remove cosmic rays and driz¬ 
zled to combine dithered exposures (Lll). Fig. [l] shows 
a time series of the combined ’R’ and ’B’ band images, 
illustrating the evolution of the supernova during the last 
20 years. 

The light curve of the ring was measured using an aper¬ 
ture in the form of an elliptical annulus, with semi-major 
axis 0.5" (1.1") for the inner (outer) boundary and an ax¬ 
ial ratio of 0.78. Corrections were applied as described 


Table 1 

New UVES and HST observations. 


Date 

EpoclF] 

(days) 


Exposure 

(<0 


UVES/VLT 


Wavelength range 


Seeing 



(nm) 


(arcsec) 

2008-11-22 - 

7,943- 

303-388+476-684 

9000 

0.8-1.0 

2009-02-08 

8,021 

373-499+660-1060 

11250 

0.8-1.3 

2010-10-18 - 

8,638- 

303-388+476-684 

9000 

0.7-1.1 

2010-11-15 

8,666 

373-499+660-1060 

11250 

0.6-1.1 

2011-10-05 - 

8,990- 

303-388+476-684 

9000 

0.8-1.1 

2011-12-04 

9,050 

373-499+660-1060 

9000 

0.6-1.4 

2013-10-18 - 

9,734- 

303-388+476-684 

8930 

1.2-1.7 

2013-12-20 

9,797 

373-499+660-1060 

8950 

0.5-1.0 

HST 


Filter 



2011-01-05 

8,717 

F438W 

1400 




F625W 

1000 


2013-02-06 

9,479 

F438W 

1200 




F625W 

1200 


2014-06-15 

9,973 

F438W 

1200 




F625W 

1200 


2014-06-20 

9,978 

F502N 

3080 


2014-06-21 

9,979 

F645N 

2880 




F658N 

2880 



a Since 1987 Feb 23 

in Lll in order to account for differences between the 
different filters/instruments as well as the loss of charge 
transfer efficiency in the late WFPC2 images. The re¬ 
sulting light curves are shown in the top panel of Fig. [2] 

From the images we can also determine the positions 
and fluxes of the hotspots. To define the hotspots, we 
used an ACS image taken in the F625W filter on Decem¬ 
ber 6, 2006, which provides the best spatial resolution 
close to the time when the ring was brightest and when 
all hot spots are present (Fig. [Tb . A two-dimensional 
Gaussian was fitted around eacnlocal maximum in the 
ring, resulting in the identification of 28 hotspots (see 
Fig. §. The distance of each hotspot from the center 
of the ring was measured by fitting a one-dimensional 
Gaussian to the intensity along a radial ray at the angu¬ 
lar location of the spot in question. In images taken prior 
to the onset of a hotspot, this method returns the radial 
location of the ring. Fig. [3] shows the time-evolution 
of the spot positions, given as a distance to the mean 
position of the ring before the appearance of the spots. 

The velocities for each spot were determined from a lin¬ 
ear fit to the positions after the impact, seen as a sudden 
decrease in the radius. For some of the spots there are 
indications of a decreasing velocity after ^ 8, 000 days. 
Given the errors in the measurements we have, however, 
limited the analysis to a linear least square fit from the 
time of the impact to the last observation. The resulting 
velocities range between 180 kms -1 and 950 kms -1 with 
a mean at ^540 kms -1 . We estimate the errors in the 
velocities to be ^ 50 kms -1 . 

The flux from each spot was determined from the 
FWHM and peak of the Gaussian. The result for each 
epoch is shown in the middle panel in Fig 3j corrected 
for differences in filters and instruments in the same way 
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Figure 1. Evolution of the ring collision from 1994 to 2014 (days 2,270 - 9,975) from a combination of HST B- and R-band images. The 
brightness of the ring has been reduced by a factor of 20 by applying a mask to the images. This makes it possible to see the morphology 
of the ring at the same time as the faint ejecta and regions outside. The new spots are better seen in the difference images in Fig. 4. The 
flux scale is otherwise the same for all panels. The size of the field is 2.1" x 1.8". North is up, east is left. 


as the total light curve. 

To complement the imaging high resolution spectra 
were obtained at regular intervals with VLT/UVES be¬ 
tween December 2000 and December 2013. These spec¬ 
tra allow us to separate the contribution from shocked 
and unshocked gas. Details of observations are given in 
Groningsson et al. (2008b) for observations up to 2007 
and later in Table [l] Reductions are similar, using the 
updated UVES pipeline. The 0.8” slit was centred at the 
supernova and oriented at a P.A. of 30°. 

We estimate the accuracy of the fluxes from the UVES 
spectra by comparing them to those from HST imaging 
in R and B bands, taking the slit width and seeing into 
account. We find that the UVES fluxes are within 20% 
of the HST fluxes, and that they also follow the same 
trend in time (Fig. J2|. The errors in the fluxes for most 
of the lines are dominated by this systematic error. How¬ 
ever, for Ha at the last epochs the separation between 
the increasingly faint narrow component from the much 
brighter shocked component becomes increasingly diffi¬ 
cult. This may introduce a systematic underestimate of 
the flux from the unshocked spots by up to ^ 50% for the 
last observations, resulting in a similar increase in these 
fluxes. The Ha light curve would then be marginally 
flatter, but the general strong decrease should not be 
affected by this. 


3. RESULTS 

The light curves of the ring in the upper panel of Fig. 
[ 2 ] are characterized by a slowly decreasing part up to day 


~ 5, 000, followed by a fast rise, which reached its peak 
at day ~ 8 , 000, after which they have decayed. The B- 
band is dominated by emission lines from H 7 , HS and [S 
II] A4069, while the R-band is dominated by Ha, [N II] 
AA6548, 6583, and [O I] AA6300, 6364. The initial de¬ 



ducing the hotspots, which brightened with time. This 
phase ended ~ 8 , 000 days after the explosion and the 
hotspots are now fading in our observations that extend 
to 9,975 days. 

The turn-on of the hot spots occurred between 5000- 
6000 days at the same time as the emission peaks moved 
inward (Fig. [3|. These dense ’protrusions’ may have 
recombined at early time due to their high density and 
only became visible when hit by the shock. The flux from 
the hotspots peaked earliest on the east side of the ring. 

The high-resolution spectroscopy from the UVES ob¬ 
servations allows us to separate the contribution to 
the light curve from the unshocked ring (em it ting nar- 
row lines with FWHM ^ 10 km s _1 ( Groningsson et al 


2008b)) and the shocked hotspots (emitting lines with 


typical velocities ~ 300 kms -1 ). In the middle panel of 
Fig. [2] we show the light curves of three different lines 
from the shocked gas (Ha, [Fe XIV] A5303 and [O III] 
A5007), which all show a similar evolution as the HST 
light curves. In the lower panel we show the fluxes of 
the strongest narrow lines from the unshocked gas (Ha, 
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Figure 2. Upper panel: R and B band light curves of the full 
ring from HST photometry. Middle panel: Light curve in individ¬ 
ual emission lines from the shocked ring from spectroscopy with 
VLT/UVES, using a 0.8 // slit. Note the different scales on the ver¬ 
tical axis in the panels and the scaling factors for the B-light curve 
and the lines. Bottom panel: The same for the narrow Ha, [N II] 
A6583 and [O III] A5007 lines from the unshocked clumps. 


[N II] A6583 and [O III] A5007). Up to ^7,000 days the 
unshocked [O III] and Ha lines increase in flux, caused 
by pre-ionization by the soft X-rays. After this epoch all 
three lines, however, decline more steeply than the lines 
from the shocked gas. 

From the UVES observations of Ha we also find that 
the shocks are radiative up to at least ^ 700 kms -1 , and 
possibly up to ^ 1,000 kms -1 , depending on the angle 


of t he shocks relative to the line of sight (Groningsson 


et al. 2008a Migo tto et ah, 2015, in prep.). Using the 


cooling time from Groningsson et al. (2008a) this indi¬ 
cates a density of up to ^ 6 x 10 4 cm -3 in the clumps. 
This can be compared to the velocities derived from the 
proper motions in Sect. [2j where the average was found 
to be ^ 540 kms -1 . This directly reflects the propaga¬ 
tion of the shocks through the clumps and and concomi¬ 
tant acceleration of the post-shock gas. Although the 
radial and proper motion velocities are in perpendicu¬ 
lar directions, the 43° inclination of the ring results in a 
similar correction, assuming that the expansion is radial 
from the explosion centre. These two methods therefore 
result in similar shock velocities. 

The most interesting recent changes are shown in Fig. 
0J with the emergence of new, faint spots, as well as dif¬ 
fuse emission, outside the ring. These new spots are most 
clearly seen in the final image from 2014, but with hind¬ 
sight, some of these are faintly visible in the 2013 image. 
The new hotspots are a factor of 10-20 fainter than the 
old hotspots. The new spots, as well as diffuse emission 
outside the ring, can also be seen in the narrow-band im¬ 
ages in Fig. [4] It is interesting that most of the new spots 



3000 5000 7000 9000 3000 5000 7000 9000 

Days since explosion 


Figure 3. Upper panel: Contour plot of the December 6, 2006 
image. Boxed regions indicate where 2D Gaussian fits were ap¬ 
plied, and the corresponding hotspot locations are identified with 
black dots. Middle panels: Evolution of the fluxes from the indi¬ 
vidual hotspots in the R band for the north-east (NE), north-west 
(NW), south-east (SE) and south-west (SW) quadrants of the ring 
numbered as in the image above. Lower panels: Evolution of radii 
from the centre (in arc secs) of the different spots as a function of 
time for the different quadrants. The radial positions are calcu¬ 
lated with respect to the average positions of the ring before the 
impact of the shock, seen as a drop at 5,000-6,000 days in the ra¬ 
dial positions. The dashed lines give a linear least square fit to the 
data after the impacts of the different spots used for the velocity 
determinations. 
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Figure 4. Upper panel: HST images of SN 1987A with WFC3 in the F625W filter on days 8,717 (2011-01-05), 9,480 (2013-02-06), and 
9,975 (2014-06-15). Middle panels: Difference images of 2013-2011, 2014-2013 and 2014-2011, respectively. Blue means fainter and red 
brighter. Inside the ring the asymmetric ejecta can be seen. Note the gradual appearance of several small spots, as well as diffuse emission, 
outside the inner ring in the south-east. New spots can also be seen in the north-east, as well as to the north-west. Note that the two 
new spots in the north-west do not coincide with the stars labelled in the upper left image. The residuals in the south-east corner of the 
difference images are due to a saturated star and the radial streak north-west of the ring in the two first difference panels is due to a 
diffraction spike of a star outside the field shown. Lower panels: Narrow band images in the F502N ([O III]), F658N ([N II]) and F645N 
(continuum) filters. The north-south asymmetries in these images are due to lines being shifted in/out of the narrow band filters.The field 
size is 3.0" x 2.4". 

are in the same region as where the faintest hotspots in 
the ring are present. The diffuse emission is likely to be a 
combi nati on of Ha from the reverse shock region (France 
et al.||2015) and gas producing line emission outside the 
ring. The latter may be a result of photoionization of 
the wind from the red supergiant stage of the progenitor 
by X-rays from the shock. 


This may explain the still increasing s oft X-ray emission 
( Maggi et al.|2012 Helder et al. | 2Q13 ), in contrast to the 


optical. Dewey et at 


. (2012|r 

illingTac 


4. DISCUSSION 

To interpret these observations we consider the ram 
pressure acting on the clumps, given by P oc p e jecta^ 2 
y-n+2^-3 ^ ^n—5^ w p ere n i s the power law index of the 
ejecta density profile, p e jecta H _n , V is the ejecta ve¬ 
locity at the reverse shock, assumed to be at a constant 
radius, and t is the time since explosion. The pressure 
will increase with time as long as the reverse shock is in 
the steep part of the ejecta density profile. This pres¬ 
sure drives transmitted shocks into clumps of gas in the 
circumstellar ring. 

If the density of a clump is sufficiently high, radiative 
cooling will cause the clump to collapse to an even higher 
de nsity . The resu lt will be an optically emitting hotspot 
(Pun et al. 2002). If instead, the density of a clump 
falls below the threshold density for radiative cooling, 
the gas behind the transmitted shock will radiate X-rays 
and very little optical radiation. This threshold density is 
a sensitive function of the velocity, VAump, of the shock 
entering the clump. Thus, at any given time, the re¬ 
quirement for radiative cooling acts as a high contrast fil¬ 
ter, limiting optical radiation to only the denser clumps. 


find from modeling of the 
X-ray observations” filling factor of ^ 30% for the dense 
clumps. 

If the radiative shocks can be characterized by one 
velocity the luminosity will depend on time as L ^ 

/JclumpVuimp- 4 « Pc\lmp ti3n ~ 15)/2A ^ A where 

Pciump is the density and A is the cross sectio n of the 
clumps and w e have assumed n « 8.6 (Shigeyama & 
Nomoto||1990|). Explosion models show that the reverse 
shock i s still likely to be in the steep density part of the 


ejecta (Dewey et al. 2012 Fransson et al. 2013). Fur¬ 


ther, the velocity for which the shocks are radiative is 
not likely to decrease because the density is increasing 
as long as n > 3. The most likely reason for a decrease 
in the luminosity is therefore that the area of the radia¬ 
tive shocks is decreasing. Alternatively, or in addition, 
as the driving pressure increases, the transmitted shock 
velocity may exceed the threshold for radiative cooling. 
Likewise, the threshold density for a clump to accommo¬ 
date a radiative shock will increase with time. In either 
case it will lead to the clumps getting dissolved. 

From the fading and acceleration of the hotspots we 
conclude that the clumps are destroyed by the blast wav e, 
in agreement with hydrodyn amic simulations (Borkowski 
Pun et ar||2002 ). The fast decay of the nar 


et al. 


1997 


row lines from the unshocked gas can be explained as a 
result of the non-radiative shocks, which traverse lower 
density gas, replacing the narrow line emission with soft 
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X-ray emission, in combination with a decreasing emis¬ 
sion measure of the pre-shock gas in the clumps. 

The range of densities in the unshocked ring has been 


estimated to 1 x 10 2 cm 3 up to 3 x 10 4 cm 3 (Mat- 


tila et al. 2010), while the s hock speed from the X-ray 
imaging is ^ 1,820 kms -1 (Maggi et al. 2012 
e tal.]2 013], coming from shocks propagating m t 


Helder 
tie low- 


density component of the gas. Assuming t hat the inter¬ 


acti on with the ring started on day 5,600 (Helder et al. 
2013), the blast wave has expanded by ~ 7 x 10 itj cm up 
to day 9,975, or ~ 10% of the radius of the ring. This 
is of the same order as the distance of the new hotspots 
outside the ring. It is therefore conceivable that the new 
hotspots are a result of either the interaction of the blast 
wave with very dense clumps embedded in gas with den¬ 
sity 10 3 cm 3 outside the ring, and close to the equa¬ 
torial plane, or pre-ionization by the X-rays. The fact 
that the main part of the ring has faded most in the 
south-east, where the majority of the new spots are seen 
supports the former interpretation. 

The time scale for the complete destruction of the ring 


in the clumps ( 

Pun et al J 2002 

). The 

mass estimate for 

the ring, 5.8 x 10 2 Mq (Matt! 

a et al. 

2010 ), only refers 


mass of the hotspots and gas outside of the ring could 
be considerably larger. However, extrapolating the light 
curves in Fig. [2] shows that the ring should fade away 
between ^ 2020 (based on the individual lines) and ~ 
2030 (based on the HST light curves). The former esti¬ 
mate is likely to be more reliable, as the lines isolate the 
emission from the shocked ring component. The hotspots 
will be gradually destroyed by instab ilities a nd conduc¬ 
tion by the hot surrounding gas ( Borkowski et al.| l997 


Pun et al. 2002). The shocked gas which has cooled anc. 


increased its density by a factor 100 — 1000 may fragment 
and survive as overdense ’nodules’ in the h o t gas, unless 
heat conduction destroys them jPol udnenko et al.f 2004 
Raga et al. 2007 Ibilvia et al.||2010|). This depends sen- 


Raga et al. 2007 Silvia et al. 20101). This depends - 
sitively on the magnetic held and its orientation. Heat 
conduction or X-ray irradiation may then provide some 
energy input to the surviving fragments, but at a low 
level. The cloud crushing time is ~ 2r c i ump /U c i ump (e.g., 


Poludnenko et al. 2004). With a decay time scale of 
2000 days and shock velocity of ~ 500 kms -1 this 
corresponds to a clump diameter of ~ 0.9 x 10 16 cm. For 
the individual spots in Sect. [2] we find a FWHM of 1.5 
pixels, consistent with being unresolved with the ACS. 
This corresponds to a diameter < 0.04", or < 3 x 10 16 
cm. We do not find any evidence of a change in the 
lateral direction within ±2°. 

The structure outside th e inner ring is complex with a 
possible hou r-glass shape (Sug erman et al.| [2005 


et al. 
the s 


France 


2015), which will now be possible to probe! As 
rock progresses beyond the circumstellar ring, it 


will trace the history of mass loss from the supernova’s 
progenitor, revealing the distribution of gas that is now 
unseen, and providing useful information to discriminate 
among different models for the progenitor of SN 1987A. 

The lower general density will result in a mainly adia¬ 
batic shock wave, except for clumps with densities com¬ 
parable to the inner ring. However, the new hotspots 


we found probably only represent a small fraction of the 
mass outside the ring. We expect SN 1987A will be¬ 
come more thoroughly X-ray dominated as the youngest 
supernova remnant evolves. 
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